Local and global instabilities of flow in a flexible-walled channel

We consider laminar high-Reynolds-number flow through a long finite-length planar channel, where a segment of one wall is replaced by a massless membrane held under longitudinal tension. The flow is driven by a fixed pressure difference across the channel and is described using an integral form of the unsteady boundary-layer equations. The basic flow state, for which the channel has uniform width, exhibits static and oscillatory global instabilities, having distinct modal forms. In contrast, the corresponding local problem (neglecting boundary conditions associated with the rigid parts of the system) is found to be convectively, but not absolutely, unstable to small-amplitude disturbances in the absence of wall damping. We show how amplification of the primary global oscillatory instability can arise entirely from wave reflections with the rigid parts of the system, involving interacting travelling wave flutter and static-divergence modes that are convectively stable; alteration of the mean flow by oscillations makes the onset of this primary instability subcritical. We also show how distinct mechanisms of energy transfer differentiate the primary global mode from other modes of oscillatory instability

A two-fluid model for tissue growth within\ud a dynamic flow environment

We study the growth of a tissue construct in a perfusion bioreactor, focussing on its response to the mechanical environment. The bioreactor system is modelled as a two-dimensional channel containing a tissue construct through which a flow of culture medium is driven. We employ a multiphase formulation of the type presented by G. Lemon, J. King, H. Byrne, O. Jensen and K. Shakesheff in their study (Multiphase modelling of tissue growth using the theory of mixtures. J. Math. Biol. 52(2), 2006, 571–594) restricted to two interacting fluid phases, representing a cell population (and attendant extracellular matrix) and a culture medium, and employ the simplifying limit of large interphase viscous drag after S. Franks in her study (Mathematical Modelling of Tumour Growth and Stability. Ph.D. Thesis, University of Nottingham, UK, 2002) and S. Franks and J. King in their study (Interactions between a uniformly proliferating tumour and its surrounding: Uniform material properties. Math. Med. Biol. 20, 2003, 47–89).\ud \ud The novel aspects of this study are: (i) the investigation of the effect of an imposed flow on the growth of the tissue construct, and (ii) the inclusion of a mechanotransduction mechanism regulating the response of the cells to the local mechanical environment. Specifically, we consider the response of the cells to their local density and the culture medium pressure. As such, this study forms the first step towards a general multiphase formulation that incorporates the effect of mechanotransduction on the growth and morphology of a tissue construct. The model is analysed using analytic and numerical techniques, the results of which illustrate the potential use of the model to predict the dominant regulatory stimuli in a cell population

Modelling crystal aggregation and deposition\ud in the catheterised lower urinary tract

Urethral catheters often become encrusted with crystals of magnesium struvite and calcium phosphate. The encrustation can block the catheter, which can cause urine retention in the bladder and reflux into the kidneys. We develop a mathematical model to investigate crystal deposition on the catheter surface, modelling the bladder as a reservoir of fluid and the urethral catheter as a rigid channel. At a constant rate, fluid containing crystal particles of unit size enters the reservoir, and flows from the reservoir through the channel and out of the system. The crystal particles aggregate, which we model using Becker–Döring coagulation theory, and are advected through the channel, where they continue to aggregate and are deposited on the channel’s walls. Inhibitor particles also enter the reservoir, and can bind to the crystals, preventing further aggregation and deposition. The crystal concentrations are spatially homogeneous in the reservoir, whereas the channel concentrations vary spatially as a result of advection, diffusion and deposition. We investigate the effect of inhibitor particles on the amount of deposition. For all parameter values, we find that crystals deposit along the full length of the channel, with maximum deposition close to the channel’s entrance

Nonuniqueness in a minimal model for cell motility

Two–phase flow models have been used previously to model cell motility, however these have rapidly become very complicated, including many physical processes, and are opaque. Here we demonstrate that even the simplest one–dimensional, two–phase, poroviscous, reactive flow model displays a number of behaviours relevant to cell crawling. We present stability analyses that show that an asymmetric perturbation is required to cause a spatially uniform, stationary strip of cytoplasm to move, which is relevant to cell polarization. Our numerical simulations identify qualitatively distinct families of travelling–wave solution that co–exist at certain parameter values. Within each family, the crawling speed of the strip has a bell–shaped dependence on the adhesion strength. The model captures the experimentally observed behaviour that cells crawl quickest at intermediate adhesion strengths, when the substrate is neither too sticky nor too slippy

A note on heat and mass transfer from a sphere in Stokes\ud flow at low Péclet number

We consider the low Péclet number, Pe ≪ 1, asymptotic solution for steady-state heat and mass transfer from a sphere immersed in Stokes flow with a Robin boundary condition on its surface, representing Newton cooling or a first-order chemical reaction. The application of van Dyke’s rule up to terms of O(Pe3) shows that the O(Pe3 log Pe) terms in the expression for the average Nusselt/Sherwood number are double those previously derived in the literature. Inclusion of the O(Pe3) terms is shown to increase significantly the range of validity of the expansion

Is the Donnan effect sufficient to explain swelling in brain tissue slices?

Brain tissue swelling is a dangerous consequence of traumatic injury and is associated with raised intracranial pressure and restricted blood flow. We consider the mechanical effects that drive swelling of brain tissue slices in an ionic solution bath, motivated by recent experimental results that showed that the volume change of tissue slices depends on the ionic concentration of the bathing solution. This result was attributed to the presence of large charged molecules that induce ion concentration gradients to ensure electroneutrality (the Donnan effect), leading to osmotic pressures and water accumulation. We use a mathematical triphasic model for soft tissue to characterize the underlying processes that could lead to the volume changes observed experimentally. We suggest that swelling is caused by an osmotic pressure increase driven by both non-permeating solutes released by necrotic cells, in addition to the Donnan effect. Both effects are necessary to explain the dependence of the tissue slice volume on the ionic bath concentration that was observed experimentally

Wrinkling, creasing, and folding in fiber-reinforced soft tissues

Many biological tissues develop elaborate folds during growth and development. The onset of this folding is often understood in relation to the creasing and wrinkling of a thin elastic layer that grows whilst attached to a large elastic foundation. In reality, many biological tissues are reinforced by fibres and so are intrinsically anisotropic. However, the correlation between the fiber directions and the pattern formed during growth is not well understood. Here, we consider the stability of a two-layer tissue composed of a thin hyperelastic strip adhered to an elastic half-space in which are embedded elastic fibers. The combined object is subject to a uniform compression and, at a critical value of this compression, buckles out of the plane — it wrinkles. We characterize the wrinkle wavelength at onset as a function of the fiber orientation both computationally and analytically and show that the onset of surface instability can be either promoted or inhibited as the fiber stiffness increases, depending on the fibre angle. However, we find that the structure of the resulting folds is approximately independent of the fiber orientation. We also explore numerically the formation of large creases in fiber-reinforced tissue in the post-buckling regime

The Thermal Structure of the Circumstellar Disk Surrounding the Classical Be Star gamma Cassiopeia

We have computed radiative equilibrium models for the gas in the circumstellar envelope surrounding the hot, classical Be star $\gamma$Cassiopeia. This calculation is performed using a code that incorporates a number of improvements over previous treatments of the disk's thermal structure by \citet{mil98} and \citet{jon04}; most importantly, heating and cooling rates are computed with atomic models for H, He, CNO, Mg, Si, Ca, & Fe and their relevant ions. Thus, for the first time, the thermal structure of a Be disk is computed for a gas with a solar chemical composition as opposed to assuming a pure hydrogen envelope. We compare the predicted average disk temperature, the total energy loss in H$\alpha$, and the near-IR excess with observations and find that all can be accounted for by a disk that is in vertical hydrostatic equilibrium with a density in the equatorial plane of $\rho(R)\approx 3$ to $5\cdot 10^{-11} (R/R_*)^{-2.5} \rm g cm^{-3}$. We also discuss the changes in the disk's thermal structure that result from the additional heating and cooling processes available to a gas with a solar chemical composition over those available to a pure hydrogen plasma.Comment: 11 pages, 8 figures high resolution figures available at http://inverse.astro.uwo.ca/sig_jon07.htm

Theory of high-energy emission from the pulsar/Be-star system PSR 1259−-63 I: radiation mechanisms and interaction geometry

We study the physical processes of the PSR B1259-63 system containing a 47 ms pulsar orbiting around a Be star in a highly eccentric orbit. Motivated by the results of a multiwavelength campaign during the January 1994 periastron passage of PSR B1259-63, we discuss several issues regarding the mechanism of high-energy emission. Unpulsed power law emission from the this system was detected near periastron in the energy range 1-200 keV. We find that the observed high energy emission from the PSR B1259-63 system is not compatible with accretion or propeller-powered emission. Shock-powered high-energy emission produced by the pulsar/outflow interaction is consistent with all high energy observations. By studying the evolution of the pulsar cavity we constrain the magnitude and geometry of the mass outflow outflow of the Be star. The pulsar/outflow interaction is most likely mediated by a collisionless shock at the internal boundary of the pulsar cavity. The system shows all the characteristics of a {\it binary plerion} being {\it diffuse} and {\it compact} near apastron and periastron, respectively. The PSR B1259-63 cavity is subject to different radiative regimes depending on whether synchrotron or inverse Compton (IC) cooling dominates the radiation of electron/positron pairs advected away from the inner boundary of the pulsar cavity. The highly non-thermal nature of the observed X-ray/gamma-ray emission near periastron establishes the existence of an efficient particle acceleration mechanism within a timescale shown to be less than $\sim 10^2-10^3$ s. A synchrotron/IC model of emission of e\pm-pairs accelerated at the inner shock front of the pulsar cavity and adiabatically expanding in the MHD flow provides an excellent explanation of the observed time variableX-ray flux and spectrum from the PSRComment: 68 pages, accepted for publication in the Astrophys. J. on Aug. 26, 199